Many-body electronic structure and Kondo properties of cobalt-porphyrin molecules. 
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We use a combination of first principles many-body methods and the numerical renormalization- 
group technique to study the Kondo regime of cobalt-porphyrin compounds adsorbed on a Cu(lll) 
surface. We find the Kondo temperature to be highly sensitive to both molecule charging and 
distance to the surface, which can explain the variations observed in recent scanning tunneling 
spectroscopy measurements. We discuss the importance of many-body effects in the molecular 
electronic structure controlling this phenomenon and suggest scenarios where enhanced temperatures 
can be achieved in experiments. 
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I. INTRODUCTION 

Predicting the electronic properties of magnetic nanos- 
tructures is a challenging task as magnetism often man- 
ifests itself as a delicate interplay between spin order- 
ing surface chemistry and finite size effects. 1,2 Nanoscale 
magnetism is at the heart of several proposals for molec- 
ular spintronic devices, highlighting the need for a de- 
tailed understanding of spin-related phenomena in these 
systems. One important and especially subtle example 
of this behavior is the Kondo effect, arising from many- 
body spin interactions between a magnetic impurity and 
electrons in a metallic system. 3 

At low temperatures, the Kondo screening of the im- 
purity's magnetic moment creates a sharp resonance at 
the Fermi energy in the local density of states at/near 
the magnetic impurity, producing a zero-bias anoma- 
lous signal in transport experiments that probe the lo- 
cal structure. Scanning tunneling microscopy (STM) 
has proved to be an outstanding tool to probe signa- 
tures of the Kondo effect in magnetic adatoms on metal- 
lic surfaces. 4-7 More recently, Kondo-like signals have 
been reported in magnetic molecules adsorbed on metal- 
lic substrates, 8-12 stimulating significant advances on the 
theory side. 13-21 For both single magnetic atoms and 
magnetic molecules on metallic surfaces, the low-bias 
conductance measured by STM acquires a character- 
istic Fano lineshapc, from which the Kondo tempera- 
ture Tk is usually inferred from fittings to well-known 
equations. 22 ' 23 Such expressions, however, do not offer 
estimates for Tk, and rather treat it as a phenomeno- 
logical fitting parameter. In fact, as we discuss below, 
quantitative estimates of Tk are not trivial to obtain in 
these systems with current theoretical approaches. 

Here, we address this issue by employing a combi- 
nation of complementary numerical methods to calcu- 
late the Kondo temperature and the density of states in 
the Kondo regime for complex magnetic molecules using 
state-of-the-art theories and a minimum of information 
from experimental data. Our approach relies on a com- 



bination of two powerful techniques: the first-principles 
GW method, 24-26 one of the most advanced tools to 
describe the electronic structure of real materials; and 
the numerical renormalization-group (NRG), 2 widely re- 
garded as the most reliable numerical method to study 
Kondo correlations. Our main goals are (i) to identify 
the relevant microscopic model and obtain realistic mi- 
croscopic parameters from the GW molecular electronic 
structure and (ii) to describe the Kondo regime in a non- 
perturbative fashion, with no additional assumptions on 
the shape of the Kondo resonance. 

The strength of this merging lies on the complemen- 
tarity of these techniques. On the one hand, ab initio 
methods face formidable challenges in the Kondo regime 
since the strongly-correlated Kondo singlet state emerges 
at energy scales much smaller than the typical atomic or 
orbital energies and is beyond typical formulations of ex- 
change and correlation effects. Furthermore, the many- 
body Kondo state is not well captured by perturbative 
approaches, which sometimes yield infrared divergencies 
at energy scales of order of /cbTr-. 3 On the other hand, 
NRG, which fully describes the Kondo regime, relies on 
quantum impurity models to capture the microscopic 
low-energy features of the full system. In the case of 
organic complexes containing magnetic atoms, the pres- 
ence of molecular states, additional degrees of freedom 
and charging effects further complicates the task of iden- 
tifying the relevant low-energy microscopic model lead- 
ing to the Kondo effect. In fact, the comparatively high 
Kondo temperatures (> 100K) reported in STM mea- 
surements involving porphyrin compounds with a mag- 
netic atom 8-12 indicate that extra, molecule- and surface- 
specific features play an active role in modifying the 
Kondo state. In all, such complexity has hindered at- 
tempts to predict Kondo temperature values from first- 
principle calculations, although recent work in nanocon- 
tacts is most promising. 28 

The presentation is organized as follows: in Sec. II we 
present GW results for the electronic structure of the 
molecule. Based on these results, we propose a realistic 
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FIG. 1: Energy of electronic orbitals of TBrPP-Co in gas 
phase, obtained from first-principles GW calculations. The 
calculated Coulomb energy cost U for double occupation of 
the d 3z 2_ 1 orbital is also indicated. Orbitals in the mid- 
dle and right columns have projection on Co 3d atomic or- 
bitals indicated. The " HOMO" and "LUMO" labels refer, 
respectively, to occupied and unoccupied orbitals in the non- 
magnetic TBrPP. The Fermi energy Ef indicated is an esti- 
mate of the actual Fermi energy of adsorbed TBrPP-Co on 
copper (see text). 



microscopic model in Sec. Ill and NRG results for the 
Kondo temperature and density of states are presented 
in Sec. IV. We summarize our findings in Sec. V. 



II. MANY-BODY ELECTRONIC STRUCTURE 

In the present work, we study TBrPP-Co [5, 10, 15, 
20-Tetrakis-(4-bromophenyl)-porphyrin-Co, stoichiome- 
try CoBr4N4C44H24], which is composed by a porphyrin 
ring with four bromophcnyl groups at the end parts and 
a single cobalt atom at its center. 9,12 ' 29 The free-standing 
TBrPP-Co has S4 symmetry. The anisotropic environ- 
ment in the central cation creates a four-fold molecular 
field which splits the cation's 3c? orbitals into three nondc- 
generate levels (with d 3z 2_ 1 , d x ?_ y ? , and d xy symmetries) 
and a degenerate pair (d xz ,d yz ). The molecule adsorbed 
on a Cu(lll) surface undergoes a planar-saddle deforma- 
tion that reduces symmetry even more and breaks the 
remaining degeneracy. 12 In the following, we restrict our 
analysis to the planar configuration and provide heuristic 
arguments for the observed reduction of Tk in the saddle 
conformation seen in experiments. 

In a preliminary step, we use DFT 30 to determine 
the electronic structure of TBrPP-Co in gas phase. We 
use norm-conserving pseudopotentials 31 for the interac- 
tion between valence and core electrons, including scalar 
relativistic effects. The exchange-correlation functional 
is calculated using the Perdew-Burkc-Ernzerhof (PBE) 
paramctrization. 32 Calculations are done using the PAR- 
SEC code. 33,34 We solve self-consistently the DFT equa- 



tions in real space, on a regular grid with grid spacing 
0.25 a. it., where 1 a.u. = 0.529 A. Electron wavefunctions 
are set to zero outside a domain enclosing the molecule. 
The distance between the domain boundary and any 
atom is at least 4 a.u., thus removing any unwanted 
effects arising from the scattering of electrons off the 
boundary. The DFT results arc compatible with previous 
calculations on TBrPP-Co 35 and similar molecules. 36 

In order to obtain accurate energies of electronic or- 
bitals, we use the GW method. 24-26 This approach is 
advantageous because it provides very accurate relative 
positions of localized d orbitals and of more extended 
sp orbitals 25,37 . More specifically, most existing approx- 
imations used in DFT underestimate the value of the 
Coulomb gap, while the GW method typically yields a 
better value and predicts molecular orbitals with energy 
ordering in better agreement with experiments. 

The GW method is essentially a perturbative approach 
based on calculating the electron self-energy by summing 
only Feynman diagrams with low order in the screened 
Coulomb interaction W. At lowest order, the self-energy 
is simply £ = iGoW, which explains the name of the 
approximation. Go here is the Green's function of an ef- 
fective non-interacting system, such as the Kohn-Sham 
(DFT) system. At the next level of approximation, 
vertex corrections are added to the previous equation: 
S = IGqWY. We use a simplified form of the vertex, built 
from the DFT local density approximation. More details 
about the construction of the LDA vertex and its impact 
on the electronic structure of confined nanostructures can 
be found in Refs. 26,38 and references therein. After the 
sclf-encrgy is calculated, we diagonalize the quasi-particlc 
eigenvalue equation, 

w 

where V e is the sum of electron-ion potential (replaced 
with a pseudopotential in our calculations) and electro- 
static potential created by electrons. Self-consistency be- 
tween the Green's function and self-energy can be added 
by replacing Go with a Green's function of the real 



system: 



.37,39,40 



iGWT, and repeating the calcula- 



tions of self-energy, polarizability, and Green's function 
until self-consistency among these quantities is found. 
We do not pursue full self-consistency because it is not 
clear how the vertex function should be updated dur- 
ing the process. Instead, we perform only one self- 
consistency cycle and assume that self-consistency is usu- 
ally obtained after a small number of iterations, if the 
starting electronic structure is similar enough to the con- 
verged electronic structure. 37,41 The results presented in 
Figs. 2 and 3 were obtained using the fa orbitals (see Eq. 
(I)- 

Owing to the simplicity of the subsequent Anderson 
model and to the fact that spin transport is dominated 
by interactions in the vicinity of the central atom, we 
carry out the GW calculations assuming that electronic 
screening in TBrPP-Co originates predominantly from 



3 



the surrounding porphyrin molecule, with very small con- 
tribution from the central cobalt atom. Following this 
assumption, we compute the screened Coulomb interac- 
tion W of the non-magnetic molecule TBrPP-Ca (with 
calcium in the center of the porphyrin ring) and use it in 
the calculation of TBrPP-Co. 

Several features of TBrPP-Co justify this methodol- 
ogy: (i) Cobalt 3d orbitals are localized and polarize 
at an energy much higher than the other valence elec- 
trons in the molecule. This is a direct consequence of 
spatial localization of those orbitals and of small hy- 
bridization between Co 3d orbitals and orbitals of the 
TBrPP compound, (ii) Co 3d orbitals retain most of 
their identity (i.e., have a high projection on atomic Co 
3d orbitals) even when surrounded by the non-magnetic 
TBrPP structure. From that, we infer that chemical 
bonds between nitrogen atoms and the central atom are 
weak, (iii) The weight of 3d orbitals in the overall screen- 
ing of TBrPP-Co is small as TBrPP has many more 
valence electrons than Co. Moreover, the properties of 
those electrons (energy, spatial distribution, electric sus- 
ceptibility) are very robust with respect to the type of 
metallic atom used, as we have checked with DFT-GGA 
calculations for compounds of the series TBrPP-M (M 
= Mn, Fc, Co, Cu, Ca), which show similar electronic 
structure and also very similar static susceptibility. For 
TBrPP-Ca, we computed the screened Coulomb inter- 
action within the time-dependent local density approxi- 
mation (TDLDA), 26 summing 1875 virtual orbitals. The 
energy difference between the highest virtual orbital and 
the LUMO is approximately 28 eV. 

As Shirley and Martin reported, 42 3s and 3p orbitals 
correlate strongly with 3d orbitals. Thus, we built a pseu- 
dopotential for Co that retains the orbitals 3s and 3p in 
the valence shell 42 and used this configuration in all GW 
calculations. In the calculations of TBrPP-Ca, we use a 
grid spacing of 0.35 a.u. 

Fig. 1 depicts the relative position of electronic or- 
bitals of gas phase TBrPP-Co, obtained from our first- 
principles GW calculations. Several orbitals around the 
highest occupied molecular orbital (HOMO) are strongly 
localized on the cobalt site and are easily identified with 
atomic 3d orbitals. All of them are occupied except for 
the d x 2_ y 2, which is completely empty, and the d 3z 2_ 1 , 
which is populated by one unpaired electron. The z direc- 
tion is perpendicular to the plane of the molecule. The 
intrinsically molecular orbitals (i.e., orbitals with weak 
3d character) are non-magnetic. The partially occupied 
orbital d 3z 2_ 1 gives a net spin of one half, J 2 = 3h 2 /A, 
to the molecule. This orbital, responsible for spin mag- 
netization, is oriented along the direction perpendicular 
to the plane of the porphyrin ring, giving a strong cou- 
pling with both the surface and the STM tip. These two 
features confirm that the d 3z 2_ 1 orbital is the origin of 
spin-dependent transport in STM experiments. 

We should emphasize the distinctions between the GW 
results presented in Fig. 1 and the corresponding calcu- 
lation within DFT-GGA (not shown). The most signifi- 



cant difference between the predictions of DFT-GGA and 
GW results are in the energy difference between occupied 
molecular orbitals and unoccupied molecular orbitals. 
This is a manifestation of the "band-gap underestima- 
tion" , which is systematically observed when one com- 
pares these two theories. In TBrPP-Co, DFT-GGA pre- 
dicts the occupied orbitals to have energies much higher 
than the GW estimate. As consequence, the HOMO- 
LUMO gap changes from 2.0 cV (GW) to 0.55 eV (DFT- 
GGA). Additionally, the value of the majority- minority 
splitting in the d 3z 2_ 1 orbital is U = 3.00 cV within 
DFT-GGA and U = 6.30 eV within GW (Fig. I). 43 

Our calculations indicate that, although the d 3z 2_ 1 
orbital is localized around the central atom, it would 
actually penetrate deep into the surface upon deposi- 
tion (sec Fig. 2). This is a consequence of the strong 
anisotropy of the d orbital and its orientation perpendic- 
ular to the plane of the molecule. Some non-magnetic 
orbitals around the Fermi energy are also extended along 
the z direction. In fact, Fig. 2(a) shows that the HOMO- 
1 (second highest occupied orbital), HOMO and LUMO 
extend in the out-of-plane direction over a similar range 
as the d 3z 2_ 1 orbital, indicating a similar penetration 
depth into the surface. 

We quantify the extension of each orbital in the z direc- 
tion by computing its accumulated density at a distance 
D from a plane containing the center of the molecule: 

h{D) = [ d 2 r /°°dz|^(r,z)| 2 (2) 

J A J D 

where i labels the orbital. GW results for Ii(D) for the 
HOMO-1, HOMO, LUMO and the d^_ x are shown in 
Fig. 3. DFT-GGA results are also shown for comparison. 

The fact that orbitals with comparable extent and dif- 
ferent angular character exist is a strong indication that 
the high Kondo temperatures and Fano lineshapes ob- 
served in experiments 9 ' 12 are the result of scattering in- 
terference involving the polarized, localized orbital and 
non-polarized orbitals in the molecule. 

As the molecule is adsorbed on the surface, hybridiza- 
tion between molecular and surface states will generally 
lead to charge transfer from the surface to the molecule. 
To verify the effect of such transfer on the levels involved 
in the Kondo effect, we performed DFT (PBE) calcula- 
tions for the Co-porphyrin molecule on a monolayer of 
Cu(lll). In these calculations, we also used the PAR- 
SEC code but now with periodic boundary conditions 
along all three spacial directions. The monolayer was 
modeled as an hexagonal superccll with 81 Cu atoms on 
the xy plane and 26 a.u. separating one monolayer from 
its neighboring images. 

We observe that the system relaxes to a configuration 
where the molecule-surface hybridization is stronger for 
levels localized in the bromine rings, whose energies lie 
well below the Fermi level. Thus, the relative position 
of the energy levels involved in the Kondo effect (i.e., 
those close to the Fermi energy) remains essentially un- 
changed once the molecule is adsorbed. These include the 
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FIG. 2: (color online) (a) Probability distribution of GW 
molecular orbitals H0M0-1, HOMO, LUMO, and d 3z2 _ l on 
a plane perpendicular to the plane of the porphyrin ring and 
crossing neighbor atoms N-Co-N. Hot colors (red, yellow) in- 
dicate high probability. Cold colors (green, blue) indicate de- 
creasing probability. The lowest amplitude blue background 
was replaced with white, therefore the extent of the dark blue 
"halo" is somewhat arbitrary. The TBrPP-Co molecule is 
shown in (b) front view and (c) side view. Atom colors are 
coded as: blue (cobalt), green (nitrogen), red (bromine), grey 
(carbon) and white (hydrogen). 

cobalt-like molecular levels and the HOMO/LUMO pair 
depicted in Figs. 1 and 2, for which the effect of charge 
transfer is chiefly a uniform shift of the Fermi energy 
of the system (as the chemical potentials of the molecule 
and the surface are aligned) , while the relative energy po- 
sitions of these molecular states and, most importantly, 
the net magnetization remains unchanged. 

We believe such a robustness of the magnetization 
arises from the relatively large separation in energy (> 
0.5 eV) of the Cobalt 3d orbitals (Fig. 1). The perturba- 
tion induced by the surface potential on the cobalt atom 
is small in this energy scale and thus the occupation of 
the orbitals does not change significantly. 

III. EFFECTIVE ANDERSON MODEL 

Using information from GW calculations, we construct 
a microscopic model with realistic parameters for the sys- 
tem and calculate the properties in the Kondo regime 
using NRG. In STM experiments with the STM tip lo- 
cated on top of the Cobalt atom, the typical width of 
the Kondo signal corresponds to small surface-tip bias 
values as compared to the characteristic spacing of elec- 
tronic levels in the molecule. In this regime, electronic 



FIG. 3: (color online) Accumulated density of orbital i (see 
Eq. 2) for orbitals d Sz 2^, HOMO-1, HOMO and LUMO of 
TBrPP-Co in gas phase. Thick lines with hlled symbols show 
the accumulated density calculated within the GW approx- 
imation. Thin lines with "+" signs show the correspond- 
ing accumulated densities calculated within DFT-GGA. The 
inset shows the ratio between integrals Id and /homo ob- 
tained both with GW (thick line) and DFT-GGA (dashed 
line). The subsequent NRG calculations were done using the 
GW /j/%OMO ratio. 



transport is dominated by elastic processes at energies 
close to the Fermi energy, Ep. This low-energy behav- 
ior can be captured by a model that includes both the 
singly-occupied d 3z 2_ 1 level as well as a molecular orbital 
M with a strong projection along the z direction at the 
central position, and with energy closest to Ep (i.e., the 
HOMO in Fig. I). 44 

The Hamiltonian for such an Anderson-like model is 
written as: 

H = H lno i + if sur f + -^coupling ■ (3) 

The first term describes the molecular levels involved in 
transport. For the description of the Kondo regime, the 
geometry of the molecule+surface+tip configuration fa- 
vors orbitals with strong z projections which are either 
singly occupied (giving the molecule an effective mag- 
netic moment) or with an energy relatively close to the 
Fermi level (contributing to electronic transport). We 
thus have 

H mo \ = ^2 E d h da + U n^hdi + E MnM<r , (4) 

a Ma 

with Ei being the GW energy of the d^ z 2 -i orbital (i = d) 
or some molecular orbital (i = M), U represents the 
Coulomb energy cost for double occupation in the d level 
and hie is the occupation number of electrons with spin 
a in orbital i. 
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The remaining terms in Eq. (3) represent the metallic 
states on the surface and the molecule-surface couplings: 

H sm { = ^ fk^-ker, 
k<r 

^coupling = ^2 (^2k4a C k<r + n - c -) + 
k 

(VjtfkC^Cka + h- c.) , (5) 

M,k 

where the operator c[ a (cio-) creates (destroys) and elec- 
tron on orbital i and Vik = ((f>i\H\ip-k) is the hybridization 
matrix element between orbital i and the Bloch state with 
vector k in the metal. 

The effects of the noninteracting molecular level terms 
(i = M) on the low-energy properties has been discussed 
in detail in Refs. 45 and 46 in the context of quantum dot 
systems, including the possibility of a "pseudogap Kondo 
effect" for Em = Ep. 45 ' 46 In the present case, we work on 
a regime where \E M — E F \ <C T K so that the effects of the 
extra level is essentially to introduce corrections to the 
Kondo temperature as well as adding an extra resonance 
to the LDOS. 

Model parameters directly obtained from first- 
principles calculations include the on-site Coulomb in- 
teraction in the d 3z 2_ 1 orbital (U = 6.3 eV) and the rel- 
ative energy separation \Ed — Em\ = 1.8 eV between the 
dj, z -2_i and the molecular orbital. The molecular-surface 
couplings, which give rise to the level linewidths Yd and 
T m , can be estimated as a function of the central atom- 
surface distance D by combining both GW results for the 
orbital spatial probability distributions and experimental 
estimates from Refs. 8,9,12, as follows. 

We consider two main contributions to the molecular 
level linewidths: an intrinsic broadening M (arising 
from electronic scattering processes within the molecule 
and obtained directly from the GW calculations) and a 
contribution 7j from the hybridization with both sur- 
face and bulk metallic states (described by -//coupling 
in Eq. (5)). The key approximation involves writing 
Ti w rp w + 7i, 47 where the surface contribution 74 is 
directly related to the hybridization matrix elements: 

7i = -Im limV lVlkl \ ■ ■ (6) 
'' k F ~ £k 111 

In general, ji will depend on the details of the overlap 
integrals VuM)k.i particularly on the distance between 
molecule and surface. While such calculations are pos- 
sible for single atoms, 13,17 they are numerically very de- 
manding for a molecule such as TBrPP-Co. 

Instead, we take the following approach: we assume 
that the relevant important metallic (bulk and surface) 
states leading to the Kondo effect have long wavelengths 
so that the surface can be replaced with an electron 
gas occupying the half-space z > D. Thus, assuming 
the metallic states to be extended and spatially homoge- 



neous, 7, on each orbital can be related to the molecule- 
surface distance D by taking 7i ex Ii(D) given by Eq. 
(2). 

We estimate the proportionality constant by choosing 
Td to be of the order of the d-level broadening mea- 
sured in low-temperature tunneling spectroscopy for pla- 
nar TBrPP-Co on Cu(lll) 9,12 (~ 0.3 eV) at a distance 
D ~ 3A, taken from Ref. 8. This sets Y d = 0.3 eV which, 
together with the calculated value for the broadening for 
the isolated molecule Tjf w = 0.15 eV, gives 7^ = 0.15 eV 
at D = 3 A. From this parameter, the ratio jd/^fM for a 
given distance D can be obtained by calculating the ratio 
between the overlap integrals Id/lM, given by Eq. (2). 

IV. KONDO PROPERTIES 

The properties of the two-orbital Anderson model 
given in Eq. (3), consisting of an interacting d-level (fi- 
nite Coulomb interaction U) and a noninteracting molec- 
ular orbital coupled to a metallic channel, can be ob- 
tained with NRG. For the NRG calculations, we write the 
Hamiltonian given in Eq. (3) as a single-channel Ander- 
son impurity model with an energy-dependent hybridiza- 
tion function 45,46 that includes the information on the 
molecular orbitals. 

We used a NRG discretization parameter A = 2.5 
and kept up to 1000 states (not including SU(2) spin 
degeneracy) on each iteration; we use a wide-band ap- 
proximation, which assumes <C i'band, where -Dband 
is the bandwidth of the continuum of metallic states 
(which includes contributions from both surface and 
bulk states). To check this assumption, we performed 
DFT calculations for bulk Cu, which give £>band ~ 10 
eV. The molecular density of states (DOS), given by 
(— Tr~ 1 )Im[Gd{ijj) + Gm was obtained from the inter- 
acting Green's functions G d ( M )(u) = ((c rf(M) : 

calculated from NRG spectra. 45,46 The DOS calculated 
in such manner contains spectral information from both 
"d" and "M" levels, which, based on their wave function 
z-extension (Fig. 3), should yield comparable contribu- 
tions to the STM signal. 

In the calculations, the orbital energies Ed and Em en- 
tering the model are referenced to the Fermi energy Ep 
at the metallic surface. While DFT gives a crude esti- 
mate for Ep (Fig. 1), we should stress it is numerically 
challenging to extract a more reliable estimate for Ep 
from a corresponding GW calculation. Instead, we chose 
to treat Ed — Ep as a free parameter. Thus, the surface- 
molecule charge transfer, identified in the DFT calcula- 
tions as the mechanism controlling the relative position 
of the molecular orbitals to the Fermi energy, can be 
studied by changing Ed and Em while keeping Ed — Em 
constant, which is equivalent to effectively varying the 
chemical potential in the molecule. 

Figure 4 shows the Kondo temperature Tr- as a func- 
tion of Ed — Ep (Fig. 4a) and of the central atom-surface 
distance D (Fig. 4b), for different Ed — Em values, as 
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FIG. 4: (color online) Kondo temperature Tk (obtained from 
NRG calculations with GW-based parameters) versus E,i — Ef 
(a) and molecule-surface distance D (b). 

shown. We find that the calculated Kondo temperatures 
are within the range of experimentally observed Kondo 
temperatures in similar systems 8 ' 9,12 for D as 3 A, and 
charging —0.8 < Ed — Ep < —0.6 eV (note that, in this 
regime, orbital M becomes the LUMO). For the value 
Ed — Ep = —0.7 eV reported in Refs. 9,12, we obtain 
Tk ~ 140 K for D « 3A, in the same range as the re- 
ported experimental results. 

Overall, we obtain an exponential dependence of Tk 
with both D and Ed — Ep. Such behavior of Tk is not 
surprising: for Ym "C Yd, the many-body Hamiltonian 
reduces to the usual single impurity Anderson model and 
one expects Tk to follow an exponential behavior: 3 

T K - ^J UTd c -T\E d -E F \\E d -E F +U\/2Ur d ^ 

The rapid decrease of Tk with D follows naturally from 
the dependence of the accumulated density Id (and hence 
Yd) on D, as depicted in Fig. 2. 

Interestingly, both exponential behaviors can, in a 
sense, compensate each other if a conformational change 
in the molecule decreases D and Ed at the same time, re- 
sulting in a weaker variation in Tk- In fact, this explana- 
tion is consistent with the behavior reported in Ref. 12, 
for which the structural change from saddle to planar 
conformation (reducing D) was accompanied by a slight 
reduction of 0.2 eV in Ed — Ep, leading to a moderate 
increase in the Kondo temperature of ATk/Tk ~ 0.3. 

The width of the Kondo resonance in the energy- 
resolved density of states (DOS) of the adsorbed molecule 
is proportional to the Kondo temperature. Experimen- 
tally, the molecule DOS can be directly probed with STM 
by suppressing the direct tunneling from tip to surface. 7 
We have calculated the DOS from the Green's functions 
for the model with NRG. Figure 5 depicts the DOS for 
different charge transfer (Em — Ep) and distance (D) 



FIG. 5: (color online) Molecular density of states (DOS) sys- 
tem for different molecular charging (a-c, D = 3 A) and 
molecule-surface distances (d-f, Em—Ef = 1.1 eV). A Kondo 
resonance at the Fermi energy and the HOMO and LUMO or- 
bitals are shown. The width variation of the Kondo resonance 
(~ Tk) is consistent with the results in Fig. 4. 



values. In all cases, the molecular LUMO resonance is 
prominent, as well as the Kondo peak at the Fermi en- 
ergy. The Ed level is present but it is not as prominent. 
The width of the Kondo peak increases substantially for 
smaller distances and changes with charging, correlating 
well with the behavior of Tk in Fig. 4. Notice that the 
Kondo feature can be extremely sharp (a couple of mil- 
livolts across) for larger molecule-surface distances (Fig. 
5d). 



V. SUMMARY 

In summary, we combine first principles GW calcula- 
tions with the numerical renormalization-group method 
and calculate Kondo temperatures of TBrPP-Co, a mag- 
netized porphyrin compound, adsorbed on copper. This 
allows for a quantitative comparison with experimental 
results, provided that two crucial parameters are known: 
the separation between molecule and surface atoms, and 
the amount of charge transfer to or from the molecule 
(or, alternatively, the energy of the HOMO relative to 
the Fermi energy of the surface). We find that our ab 
initio treatment beyond DFT is essential, as the contri- 
bution from dynamical correlations to the molecule pa- 
rameters (e.g., U , orbital energies and hybridizations) are 
important. 

Our first-principles GW calculations indicate that 
there are two kinds of orbitals in the vicinity of the 
Fermi energy with a strong overlap with the surface: 
a spin-polarized d 3z 2_ 1 orbital, which gives rise to the 
Kondo effect, and non-polarized molecular orbitals. Both 
types have strong vertical projections and should dom- 
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inate the observed STM conductance signal. Numeri- 
cal renormalization-group calculations for a microscopic 
model built based on these GW results yield Kondo tem- 
peratures in good agreement to those observed in re- 
cent STM measurements. 9 ' 12 Additionally, the calculated 
DOS shows both Kondo and molecular resonances con- 
sistent with the tunneling spectroscopy data. 

We find, not unexpectedly, that the Kondo tempera- 
ture and the DOS features depend strongly on the inter- 
play between two parameters: the position of the par- 
tially occupied spin-polarized orbital with respect to the 
Fermi energy, which increases Tk as it approaches that 
energy, and the wave function overlap between molecu- 
lar orbitals and surface states, which decreases Tk as it 
reduces (increasing D). Most interestingly, in a regime 
where these two parameters compensate each other, Tk 
may show non-exponential behavior providing a natural 
explanation to recent STM measurements of Tk ■ 

We should note the combination of numerical meth- 
ods employed here is quite general and can be used to 
study other mctallo-organic molecules. Examples are 
porphyrin molecules with different magnetic atoms in 
the center such as Fe, Cu and Mn. Although the low- 



energy model extracted from the GW calculations in 
such cases will, in general, involve multiple orbitals and 
higher spin configurations, such situations can be ac- 
commodated within multi-orbital extensions of the NRG 
method, 27 bringing the interesting prospect of a compre- 
hensive theoretical description of such strongly correlated 
molecular systems. 
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